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We treat the spin injection and extraction via a ferromagnetic metal/semiconductor Schottky 
barrier as a quantum scattering problem. This enables the theory to explain a number of phenom- 
ena involving spin-dependent current through the Schottky barrier, especially the counter-intuitive 
spin polarization direction in the semiconductor due to current extraction seen in recent experi- 
ments. A possible explanation of this phenomenon involves taking into account the spin-dependent 
inelastic scattering via the bound states in the interface region. The quantum-mechanical treatment 
of spin transport through the interface is coupled with the semiclassical description of transport 
in the adjoining media, in which we take into account the in-plane spin diffusion along the inter- 
face in the planar geometry used in experiments. The theory forms the basis of the calculation 
of spin-dependent current flow in multi-terminal systems, consisting of a semiconductor channel 
with many ferromagnetic contacts attached, in which the spin accumulation created by spin injec- 
tion/extraction can be efficiently sensed by electrical means. A three-terminal system can be used 
as a magnetic memory cell with the bit of information encoded in the magnetization of one of the 
contacts. Using five terminals we construct a reprogrammable logic gate, in which the logic inputs 
and the functionality are encoded in magnetizations of the four terminals, while the current out of 
the fifth one gives a result of the operation. 

PACS numbers: 



I. INTRODUCTION 

Electrical spin injection (transfer of spin polarization by 
electrical current) from a ferromagnet into a paramag- 
net was first achieved in junctions between metals by 
Johnson and Silsbecj^i^ The spin injection into semicon- 
ductors has proven to be a harder task.'^ In late 90s a 
successful spin injection at low temperatures from Mn- 
doped diluted ferromagnetic semiconductors^i^ii gave 
new impetus to the field of semiconductor spintronics. 
Injection from ferromagnetic metals at temperatures up 
to the room temperature followed soon afterwardsj ^i^i^° 
Initiallj*^ the injection had quite low efficiency, which 
was was later increased j^iliii^ with maximum reported^^ 
value of 30%. The rise in spin injection efficiency was 
achieved by a proper doping of the metal/semiconductor 
interface,— d^d^ Spin extraction was also seen through 
optical measurement of spin accumulation in forward- 
biased MnAs/GaAs junctionii^ 

During the last two years, there was a tremendous 
progress in both spin injection and extraction in Fe/GaAs 
structures. The spin accumulation due to both spin 
injection from Fe and spin extraction from GaAs into 
Fe (a depletion of spins which can move more easily 
into the magnet) were imaged by Kerr spectroscopyJ^ 
Soon afterwards the spin accumulation in the semicon- 
ductor near the junction with a magnet has been sensed 
electrically ) ^^d^ proving that the current through the 
metal/semiconductor junction depends on spin polariza- 
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tion of electrons inside the semiconductor. 

In these experiment o^^d^d^ an unexpected sign of spin 
accumulation near the drain was seen. For source and 
drain with parallel magnetizations, and for potential 
drops at the interfaces small compared to the energy scale 
on which the spin-projected densities of states in the fer- 
romagnets change significantly, spins of opposite direc- 
tions should accumulate near source and drain contacts. 
Contrary to this expectation, which is a consequence of 
time-reversal symmetry for elastic tunneling between the 
two materials?^ the observed spin accumulation near the 
drain had the same sign as the one near the source con- 
tact. In Sec. Ill Bi we review a theory^! of spin extraction 
which takes into account inelastic scattering through the 
bound states near the interface. These arise from the in- 
homogeneous profile of heavy n-type doping of the semi- 
conductor near the junction with the meta L^^d*^ ;^ This 
doping makes the Schottky barrier thin enough (^10 
nm) for efficient spin transport via tunneling. It also 
results in the creation of a potential well for electrons 
next to the barrier. At forward bias, the electrons from 
the bulk of the semiconductor either tunnel directly into 
the metal (elastic process), or scatter inelastically into 
the quasi-bound states in the well, and then leak out 
into the magnet. These two channels of electron trans- 
port through the interface favor opposite spin orienta- 
tions, leading to opposite signs of spin accumulation. An 
alternative theory based on first-principles calculation of 
interface electronic structure (but neglecting the bending 
of the conduction band potential in GaAs) has also been 
proposed. •^^ 

Efficient spin injection/extraction is a basic prereq- 
uisite for any kind of practical application of spintron- 
ics systems. The simplest spintronics device is a two- 
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terminal spin valve, in which a current is passed between 
two ferromagnetic contacts connected through a param- 
agnetic channel. However, in electronics three-terminal 
semiconductor devices are indispensable for their switch- 
ing (biasing the gate of a field effect transistor) or am- 
plification capabilities (driving a current into the base of 
a bipolar transistor). Many types of "spin transistors" 
have been proposed theoretically. The most famous is 
the simple "current modulator" proposed by Datta and 
Das^S, in 1990, in which the electric field of the gate 
together with Rashba spin-orbit interaction^iiS^ in the 
small-bandgap semiconductors controls the spin preces- 
sion of ballistic electrons injected and extracted by fer- 
romagnetic contacts. Despite a large experimental effort 
a conclusive demonstration of the device operation has 
remained elusive. Let us mention some of the other pro- 
posed semiconductor spin-transistors. A diffusive ver- 
sion of a Datta-Das system has been put forthiS" Mag- 
netic unipolar^^ and bipolar transistor o^^i^^'^" (both of 
which require non-degenerate magnetic semiconductors) 
have been analyzed. Another proposal was that of a spin 
transistor without any ferromagnetic elements?^ which 
relies exclusively on strong spin-orbit interaction expe- 
rienced by electrons in small bandgap materials such as 
In As. There was also an idea of bypassing the prob- 
lems with efficient spin injection and using a proxim- 
ity effect of a ferromagnetic gate j^^'^'^ Here, we review 
our theoretical work on a class of systems consisting of a 
semiconductor channel with multiple ferromagnetic con- 
tacts attached. We work in the regime of diffusive spin 
transport at room temperature, and we concentrate on 
Fe/GaAs structures. However, we do not exploit strong 
spin-orbit interaction present in III-V semiconductors; 
the spin-orbit scattering serves only as a source of spin 
relaxation. Consequently, the proposed devices are also 
suited for silicon-based systems, as the spin relaxation 
time in Si is expected to be at least an order of magni- 
tude longer than in GaAs. This makes silicon a perfect 
candidate for spintronics applications which do not rely 
on spin manipulation through spin-orbit interaction. The 
existence of spin-dependent coupling through a heavily 
doped Schottky junction between Si and a paramagnetic 
metal has been shown indirectlyj^ Recently, progress has 
been made^ in creating tunneling contacts with widely 
tunable conductance between a ferromagnet and Si, al- 
beit without showing yet spin injection. It is encouraging 
that hot electron spin injection into silicon accompanied 
by a magnetoresistive effect has been achieved recently,^^ 
as well as spin injection from iron into silicon using alu- 
minum oxide tunneling barriers i^l 

Our review is organized as follows. In Sec. |TT] we give 
a theory of spin-polarized transport through the heav- 
ily doped metal/semiconductor interface. Spin injection 
at large reverse bias is described in Sec. Ill Al For spin 
extraction we introduce in Sec. IIIBI a new mechanism 
of spin transport through the junction due to leakage of 
localized electrons from the potential well close to the 
interface, and we show that it gives an opposite sign of 



spin accumulation to the mechanism of direct tunneling 
between the bulk of the semiconductor and the metal. 
Then in Sec. Ill CI we consider a special case of an "op- 
timally doped" barrier (without a pronounced potential 
well) kept at low bias, which can be used as an electri- 
cal probe of the spin accumulation in the semiconduc- 
tor. In Sec. mil we couple the description of spin in- 
jection/extraction with the diffusive transport inside the 
semiconductor, and we explain the significance of spin ac- 
cumulation for magnetoresistance of a two-terminal sys- 
tem (a spin valve). We introduce the basic concept of 
electrical sensing of the spin splitting in the semicon- 
ductor using a ferromagnetic contact in Sec. [IVl The 
possible applications, including the Magnetic Contact 
Transistor— and a magneto-logic gate integrable into a 
large-scale circuit^ are presented in Sec. IIVI and are il- 
lustrated by theoretical calculations. 



II. SPIN TRANSPORT THROUGH THE 
SCHOTTKY BARRIER 

The Schottky barrieri^Siil between a metal and a semi- 
conductor is created by redistribution of charges in the 
space charge layer. We denote the barrier height mea- 
sured from the Fermi level of the metals by 0b, and its 
thickness as d. For uniformly n-doped semiconductor the 
barrier shape is approximately parabolic, and the deple- 
tion width d is the distance between the interface and 
the onset of the bulk flat-band region. 

For Fe/GaAs junctioni^ we will employ a value of 
0B=O.8 eV. For homogeneous doping no<10^^ cm~^ we 
have the depletion width d>lQQ nm. For such a wide bar- 
rier the tunneling current is negligible, and the current is 
due to purely classical thermionic emission*^ which de- 
pends on temperature and the barrier height, not on its 
width or shape. Even if this current is spin polarized, its 
total density is too small to create an appreciable spin ac- 
cumulation. For higher no the tunneling dominates the 
transport through the junction, but only for extremely 
high bulk doping levels (no~10^^ cm~'^) we have c?«10 
nm resulting in appreciable current densities. In order 
to achieve such thin barriers yet with the bulk of the 
semiconductor having the carrier density less than 10^^ 
cm~'^, a strongly inhomogeneous doping profile has to be 
used near the interfaccj ^'^'^^ '^ Spin injection from Fe into 
GaAs with bulk no~10"'^^ cm~'^ has been observed ^'^^i^^ 
only in such heavily doped junctions, in which the first 
15 nm of semiconductor beneath the interface is doped 
with nd=5 ■ 10^^ cm^"^ donors. 

In general, the doping of the interface results in a cre- 
ation of a potential well close to the barrier Even 
if there is no well in equilibrium, at high forward bias 
when less electrons need to be depleted from the semi- 
conductor, the well creation is inevitable. In Sec. Ill Bl 
we show that the presence of bound states in this well 
can have a profound effect on spin extraction from the 
semiconductor . 
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A. Spin injection 



Theoretical analysis of spin injection from metals into 
semiconductors has shown^^i^i^ that the junction with 
large resistance (a tunneling barrier) is necessary for 
the current to be polarized. More precisely, since the 
spin-depth conductance of the semiconductor G^c=a jL 
(with conductivity a and spin diffusion length V) is much 
smaller than its metal counterpart Gm=am/ L,n^ for spin 
injection to occur the junction conductance G has to ful- 
fill G<Gsc ■ In such a case the spin polarization of the cur- 
rent at the interface is determined by the spin-selectivity 
of the barrier, AG=G+—G-, in which Gs are the conduc- 
tances for spin s=± (along the quantization axis given by 
magnetization of the ferromagnet). This "conductivity 
mismatch" effect was actually first analyzed in 1987 by 
Johnson and Silsbee4^ From experiments on Fe/GaAs, 
AG/G<0.3 was deduced for spin injectionpii We also 
stress that although the barrier with G^Gsc gives spin- 
polarized currents, the total current density can be too 
small to create an appreciable spin accumulation in the 
semiconductor.—!^'^ A rule of a thumb is that G^Gsc 
leads to efhcient spin injection (i.e. resulting in large spin 
accumulation), but, strictly speaking, the geometry of a 
system has to be taken into account when choosing the 
optimal barrier parameters 4^ 

An important quantity in the description of spin trans- 
port is a spin-dependent electrochemical potential /^^(x). 
It is defined as 



A^sW = M^(x) - e0(x) , 



(1) 



where /i^ is the chemical potential of electrons with spin 
s=±, (j) is the electrostatic potential, and the elemen- 
tary charge e>0. The spin splitting of the electrochemi- 
cal potential, Afi—fi^ — fi^, corresponds to the presence 
of non-equilibrium spin density (spin accumulation). In 
a non-magnetic material A/x^O means An= n+— ri_7^0, 
where is the density of electrons of spin component s. 

In Fig.[T^ we show the energy diagram of the Schottky 
barrier. We define the bias eV applied to the junction as 
the difference between the average electrochemical poten- 
tial in the flat-band region /i=(/i+ -|- /Lt_)/2 and metal's 
V>0 (<0) is forward (reverse) bias corresponding 
to electrons going from (into) the semiconductor. This 
definition of V is convenient in the case of spin accumu- 
lation small enough for fig to be linearly proportional to 
the nonequilibrium parts of the spin densities Sus, see 
Sec. mil Then, because of quasi-neutrality^ we have 
(5n+ -|- Sn-=0 and fi is equal to the equilibrium chemical 
potential in the semiconductor. 

A schematic picture of spin injection and extraction 
through Schottky barriers is shown in Fig. [T]3. The junc- 
tions are much more resistive than the semiconductor 
channel, so that the electrochemical potential shows dis- 
continuities at the barriers. Because of this and the large 
difference of conductivities of the semiconductor and the 
ferromagnet, we can disregard both spatial and spin de- 
pendence of fj,s in the ferromagnet, and use a single value 
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FIG. 1: (a) Heavily doped Schottky junction at reverse bias. 
The semiconductor is non-degenerate (at room temperature) . 
fi is the average electrochemical potential at the onset of 
the flat-band region, and A/x is the spin accumulation, (b) 
Schematic picture of current flow from the ferromagnetic in- 
jector down to the ferromagnetic spin-extracting drain. The 
potential wells filled with carriers near the interfaces are cre- 
ated by the inhomogeneous doping profile (heavy doping 
near the interface). For extraction, two routes are drawn: 
a direct tunneling from the bulk of the semiconductor (elas- 
tic process) and tunneling of the electrons from the potential 
well into the magnet sustained by capture of bulk electrons 
(inelastic process). 



of chemical potential /x™. The current injection occurs 
because of tunneling of electrons from the metal into the 
semiconductor, described^fi by spin-dependent transmis- 
sion coefficient of the particle flux Ts{kx, ky), with the 
wave-vector in the direction of the interface, and k|| the 
in-plane wave-vector, which we assume conserved (spec- 
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ular transmission). 

We neglect the atomic structure of the Fe/GaAs 
interface,— 1^1^ and use a simpUfied band-structure for 
the bulk ferromagnet with a single spin-split band with 
spin-dependent Fermi velocities Vm,s- For iron, a model 
with effective mass mm equal to the free-electron mass 
and Fermi wave- vectors A:™ =1.1 ( A:™ =0.42 A~^) for 
majority (minority) electrons has been widely usedi^i^ 
Due to the assumption of specular transmission only the 
electrons having s~fm,s can tunnel into the F val- 
ley of the semiconductor's conduction band. We define 
the imaginary wave vector within the barrier n and the 
corresponding velocity Vb=hK/msc (with semiconductor 
effective mass msc)- For high barrier and msc-^nrim con- 
sidered here we have Vb^Vm,s^Vsc, with the transmitted 
electron velocity in the semiconductor Vsc- Within this 
model we obtain for the flux transmission coefficient: 

rp ^ '^m,s'^sc ^-2Kd ^ Vm.s'^sc ^-2Brf _ A -2K,d 

If we approximate the barrier by a square step of thick- 
ness d, then K=^j2msc{4>B + M ~ e.V)/h. For a triangu- 
lar barrier the expression for k has to be modified,— but 
the spin-dependent Ag factor remains the same. Conse- 
quently, electrons with larger velocity in the metal tun- 
nel more efficiently into the semiconductor. In iron this 
translates into preferential injection of majority spins, 
which is in agreement with experiments.^^ Within this 
model we also expect that if Fe is replaced by a zinc- 
blende MnAs, because of the different ratio between 
the majority and minority spin wave-vectors,^^ the spin- 
selectivity of the junction will have opposite sign to the 
Fe case. 

For large reverse bias {\eV\':S>kBT) the injected cur- 
rent does not depend on the occupation function in the 
semiconductor, since most of the electrons tunnel from 
the metal into the states at least ksT above the chemi- 
cal potential in the semiconductor (see Fig. [Dd). Up to 
a certain critical reverse bias the barrier thickness does 
not change much, only the well becomes more shallow. 
Above this critical bias the electrons start to be depleted 
from the bulk of the semiconductor. The wide depletion 
region created then was shown to be detrimental to spin 
injection.^'' Another reason for avoiding too large reverse 
biases is that hot electron injection is accompanied by 
enhanced spin relaxation in GaAsi^ 



B. Spin extraction in the presence of bound states 
near the interface 

An analogous calculation of tunneling from the 3D states 
in the bulk of the semiconductor into the metal (spin ex- 
traction) gives the same spin selectivity, so that the spins 
parallel to the minority spin in Fe should be accumulated, 
contrary to the observation.^? The experiments can be 
explained by including the presence of electrons localized 
in in the well near the interface, and considering a two- 



step process, in which tunneling of the electrons from the 
bound states into the ferromagnet is followed by vacant 
states being filled by decay of extended state electrons 
(carrier capture). 

The spin-selectivity of the junction for free and local- 
ized electrons is explained in the following way. For free 
electrons, current conservation is well defined on both 
sides of the interface region (elastic scattering). As ex- 
plained before, in reverse or in low forward bias the 'ef- 
fective velocity' in the barrier dominates and the current 
scales with the electron's velocity in the metal side. On 
the other hand, for localized electrons the conservation 
of total reflection and transmission is irrelevant. Elec- 
trons escape from the well into the vacant states in the 
ferromagnet and the transmitted current scales with the 
decay rate of the bounded wave function. In this case, 
the decay is fastest when the electron's velocity in the 
ferromagnet matches the 'effective velocity' in the well 
(being inversely proportional to the well's width). Later 
we show, that in the case of Fe/GaAs one gets an antipo- 
dal spin-behavior for free and localized carriers j^i 

We denote the bulk doping density by tt-q. The ultra- 
heavily doped region at the junction has width d and 
doping Ud- There is also a transition region of width dtr 
where the donor density interpolates between no and Ud- 
The conditions for the existence of the potential well are 
the folllowing>2i 



e^Hd 



Ud '> m , 
dtr ^ As , 



(3) 

(4) 
(5) 



where is the relative permittivity of the semiconductor 
and \b is the electron's de Broglie wavelength (typically 
~10 nm in GaAs quantum wells). The first equation 
determines the Schottky barrier thickness d, the second 
guarantees an excess of electrons in the transition region, 
and the third one is for the two-dimensional character of 
the electronic states in the well. 

In Fig. [5^ we show the results of the self-consistent 
calculation of Schrodringer and Poisson equations for 
the conduction band profile near the junction at 0.2 
V forward bias using no=3.6 • 10^^ cm~^, nd=^ ■ 10^^ 
cm^'^, rf=(itr=15 nm and €^=12. 6. These parameters 
are chosen to approximate the junctions used in the 
experiments. ^^'^^'^^ At this bias three bound states are 
present in the well. 

The process of electron escape from the well into the 
continuum in the ferromagnet is calculated in the follow- 
ing way.^^ At time i=0 the wave functions for spin s=± 
are identical, taken as the ?th bound state, and they are 
equal to zero on the metal side. The time-dependent 
Schrodringer equations is then solved numerically using 
the potential from Fig. [2^. Inside the metal (iron) we 
use the same simplified model as before. For the cal- 
culation we have used a one-dimensional box of 200 nm 
width encompassing both the well and the inside of the 
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metal, and we have used discrete transparent boundary 
conditions to prevent reflections from the edges »^ Fig.l^b 
shows the wave functions of the third quasibound state 
penetrating into the metal at times t—4Q and 400 fs. The 
electrons with minority spin in Fe have bigger penetration 
amplitude, and this behavior persists for longer times and 
for all the bound states. The escape rate is practically 
constant in time (resulting in exponential decay of the 
quasibound state) and is given by 



1 



dx\tpi^s{x,t)\^ 

semiconductor 



(6) 

Fig. [2t shows the escape rate from the first quasi-bound 
state versus the wave- vector in the metal ■ The escape 
rate has maximum when the 'effective velocity' in the well 
matches the velocity in the metal: fc„i/r7i,„^7r/(mscdtr)- 
The values of Fermi wave- vectors which we use for iron 
are on the right side of the curve, where the escape rate 
decreases with km- Such behavior for high metal wave- 
vectors agrees with the extended WKB model for the 
alpha-particle decay,^- in which the coupling between a 
quasi-bound state and the continuum scales with the in- 
verse wave-vector in the continuum. 

The spin-dependent current density J2D.S due to leak- 
age of localized electrons is proportional to the areal elec- 
tron density in the ith state hi^s (with energy above 
the metal's Fermi level) divided by the escape time: 
J2D,s'^^i fi-i,s/T^1s'^. The spin relaxation time in the well 
is of the order of tens of pSf^i whereas the escape time 
is ^1 ns, so that we have little spin accumulation in the 
well: hi^s—fii/'^- The electron which escapes from the 
well into the magnet is replenished by an electron with 
the same spin from the bulk region due to spin- conserving 
capture process of free electrons by the well, e.g. by 
emission of longitudinal optical phonons or carrier-carrier 
scattering with degenerate electrons in the well. These 
processes are much faster than the spin relaxation time in 
the welL^^ Thus, the bulk region is left with more spin-up 
(down) electrons if it provides the well with more spin- 
down (up) electrons. The necessity of the capture process 
is consistent with the longitudinal optical phonon signa- 
ture in the low temperature conductance measurement 
of Fe/GaAs by Hanbicki et alJ^ The presence of this sig- 
nature for the forward bias had been an open question, 
and our model suggests a possible explanation. 

Apart from the current J2D due to the escape from the 
well, there is always some current Jb due to direct tunnel- 
ing between the bulk of the semiconductor and the metal. 
Fig. [5Ji shows the contributions to the total current due 
to the elastic and inelastic processes versus the bulk dop- 
ing no, calculated at zero temperature. The transmission 
from the bulk states is calculated for the self-consistent 
potential using the transfer matrix method, including 
the resonant behavior of free electrons due to the well 
potential. Since this potential is weakly affected by 
changing uq as long as no<Cnd, the J2D current increases 
only by 10% in the shown range of rip. On the other hand, 
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FIG. 2: (a) The self-consistent conduction band potential in 
the semiconductor for 0.2 V forward bias with the wave func- 
tions of the three bound states, (b) Spin-dependent ampli- 
tudes of the wave function penetrating into the ferromagnet 
after 40 and 400 fs. (c) Escape rate from the third bound 
state versus the wave-vector in the metal. The upper (lower) 
mark refers to the minority (majority) electrons in Fe. (d) Ex- 
tracted current density as a function of the backround doping 
no at zero temperature. Ji, and J2D denote current due to free 
and localized electrons, respectively, and ,/[, + J2D denotes the 
sum of the two terms (the total current), (e) Spin polariza- 
tion of the total current. Panels (d) and (e) are adapted from 
Ref. 20. 



Jb depends strongly on uq: as the chemical potential in 
the bulk increases with no, the number of carriers which 
can tunnel into the metal increases. Fig. ^ shows the 
spin polarization of the total current Pj=(J+ — J-)/ J. 
The critical background doping at which Pj—Q is 1.5T0^^ 
cm"'^. This is not exactly the density at which Jb=J2D 
because \Pj^ \ and iPj^j^l are slightly different. 

Because the effective spin selectivity of the junction 
depends on the carrier density and the bending of the 
conduction band, it can be controlled by voltage in the 
electrical spin switch^ In a semiconductor channel of 
thickness h^lOO nm, a voltage applied to a back-gate 
deposited on the opposite side of the channel to the fer- 
romagnetic contact can be used to manipulate the den- 
sity of free electrons below the magnet. By varying this 
back-gate voltage Vg in a properly designed structure it 
is possible to switch between the transport regimes dom- 
inated by tunneling from the bulk and escape from the 
well, and thus change the sign of the spin accumulation 
in the semiconductor!^ 
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C. Spin-dependent conductance at low bias 

Now let us analyze the conductance of the interface with 
a very shallow well at small bias (an "optimally doped" 
barrier). We concentrate on the room-temperature case 
relevant for potential applications, and take the bulk 
of semiconductor as non-degenerate. An example of 
a doping profile yielding such an interface is the S- 
doping)^i^i^ in which a single monolayer of a semicon- 
ductor material near the interface is doped with a donor 
density impossible to achieve in the bulk material. When 
a ^-doping layer is placed at a distance do from the inter- 
face, and its planar density is n2D=€r^o{4>B + A')/47re^do5 
then the barrier shape is triangular, there is no well 
at zero bias, and the barrier width is dp. For (io=5 
nm, 0B=O.8 eV, and /x=— 0.1 eV (corresponding to bulk 
n=10^^ cm~^), we get n^D^lO^^ cm^^. 

For the "optimally doped" contact only the elastic 
transport channel is present, and we can derive an 
expression^ for the current with spin s. In the regime 
of bias \eV\, |/^™ — A'sl, and — p\ smaller than fcsT 
the formulas for js^-, fJ-s) can be linearized around 1^=0, 
and we obtain for the spin currents: 



4e2 

Gs = A^e- 



-2K,d 



no , 



(7) 

(8) 



where Gs is the barrier conductance at low bias, the spin- 
dependence of which comes from the As factor. As dis- 
cussed above, the ratio G+ / G_ is equal to the ratio of the 
velocities of carriers with different spin in the ferromag- 
net, which is approximately 2 in our model in agreement 
with spin-LED experiments.^'^"'' For no— 10^^ cm~^ and 
do~5 nm we obtain Gs of the order of 10'^ il^^cm^^. In 
the following sections, we will use the above model of 
spin-dependent properties of the junction to model how 
the spin accumulation can be sensed electrically by a con- 
tact kept at low bias. 



III. SPIN ACCUMULATION IN THE 
DIFFUSIVE SPIN VALVE 

We work in the diffusive regime, in which the spin relax- 
ation time Tsr is much longer than the momentum scat- 
tering time. Then, from the Boltzmann equation we can 
derive the spin diffusion equation^'^= '^^i^^i^^ for the non- 
equilibrium parts of the spin densities Sus- In a para- 
magnetic and non-degenerate semiconductor the electro- 
chemical potential /i^ defined in Eq. ([1]) is given by 



fig = ksT In 



no/2 + Sris 
no/2 



knT 



no/2 



(9) 



where no is the total carrier density, and the second ex- 
pression is the linear approximation valid for \ Afi\<kBT 
(equivalently \Sns\<no/2). Here we concentrate on the 



linear regime and low electric fields,— so that we can 
write the diffusion equation for the spin-splitting of the 
electrochemical potentials: 



Ms 



Ms 



fJ-s 



2L2 



(10) 



where the spin diffusion length is defined in terms of dif- 
fusion constant D and spin relaxation time by L=\/ Dtst- 
For the spin-s current in the paramagnetic semiconductor 
we have 



-V/is = (TsE + eD\/ns , 



(11) 



where the conductivity for spin s is Us—nscv with v being 
the mobility, and cr_|_~cr_ to the first approximation in 
the linear regime. The only way for the semiconductor 
to support a non-zero spin polarization of the current is 
by creating a net spin density {n^^n^ corresponding to 
A^i^O). 

In the spin valve the current is passed through a para- 
magnetic channel between two ferromagnetic contacts. 
We assume collinear magnetizations, oriented either par- 
allel (P) or antiparallel (AP) with respect to each other. 
If the distance between the contacts is smaller than spin 
diffusion length L, the spin accumulation in the channel 
depends on the alignment of the magnets. Provided that 
the spin transport mechanism is the same for injection 
and extraction, in P configuration the same spin species is 
preferentially injected and extracted. The spin accumu- 
lation has opposite signs in the neighborhood of the two 
contacts, and by lA/i'^j we denote the magnitude of spin 
splitting near the junctions. In the AP configuration, 
spins of opposite directions are more easily injected and 
extracted, resulting in large and nearly uniform spin ac- 
cumulation: |A/L("*^|3>|AyU'^|. If the mechanisms of spin 
injection and extraction are different, and the spin selec- 
tivities of injecting and extracting junctions are opposite 
as described in Sec. Ill B[ the labels P and AP refer not to 
the relative orientations of contact magnetizations, but 
to the spin accumulation patterns described above. 

The magneto-resistive (MR) coefficient of the spin 
valve is commonly defined as MR={I^—I^^)/I^, where 
jP(AP) |-j^g total current between the two terminals. A 
qualitative relation between the MR and the spin accu- 
mulation can be derived using the simple boundary con- 
ditions for currents from Eq. ([5]). 

For the channel with relevant dimensions smaller than 
L, by balancing the net spin injection into the channel 
with the spin relaxation we obtain 



AP 



MR 



AG All 
"G eV~ 



1 



V Gsc 
2AL G 



(12) 



where V is the volume of the channel and A is the area of 
the contacts, which are assumed to both have the same 
Gs and A for simplicity. The MR depends on the ratio of 
Afi^^ to the applied bias, which for small electric fields is 
just a geometry-dependent constant. For realistic param- 
eters of a Fe/GaAs spin valve we obtain M _R^(AG/G)^, 
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and a typical value of MR is about 1% after optimizing 
the system's geometry4^ In a one-dimensional channel 
of length I we have MR^{G/Gsc)-{L/l) which is quite 
small for realistic values of G and practical values of 
Z^lOO nm. On the other hand the ratio of spin splittings 
\Afj,^^ /A^^\^{2L/l)'^ can be quite large even when AIR 
is small. This large difference of spin accumulations in P 
and AP cases is not expressed in a two-terminal system 
in the most effective way by the magnetoresistive effect. 

The MR of a one-dimensional spin valve has been 
calculated analytically^^ also using the boundary con- 
ditions from Eq. ([8]). The spin valve in the lateral geom- 
etry relevant for experimentsliilS.ii^ has been analyzed 
qualitatively^2i^ and quantitatively!^ In the latter work 
an effective one-dimensional diffusion equation was de- 
rived, accurately describing the spin diffusion in a layer 
of material of thickness h smaller than the spin diffusion 
length L, and covered by contacts with junction con- 
ductances G smaller than the conductance cr/h of the 
underlying semiconductor layer. In a structure like the 
one shown in Fig. [3^, we calculate the spin diffusion by 
introducing the layer-averaged electrochemical potential 
S,s{x) in the semiconductor channel: Cs — j^ Jq dy Hs{x, y). 
By integrating out the y dependence from Eq. (|10p we 
obtain the approximate equation: 



2L' 



Of 

+ ^{is-r), (13) 

an 



where the second term on the right-hand side is present 
only under the contacts, and /i™ is the electrochemical 
potential in the ferromagnet. This equation is derived 
using the boundary condition from Eq. ([S]), and assum- 
ing small electric fields and small spin accumulations (so 
that (Ts«ct/2). For Fe/GaAs structures with Gg^lO^ 
f2~-'^cm~^ this approximate formalism gives results indis- 
tinguishable from exact numerical calculations, and all 
the results presented below are obtained using this ap- 
proach. 



IV. ELECTRICAL EXPRESSION OF SPIN 
ACCUMULATION IN MULTI-TERMINAL 
SYSTEMS 

In the previous section we have seen that in the spin valve 
the patterns of spin accumulation in the semiconductor 
are qualitatively different for P and AP configurations, 
but the MR ratio does not directly reflect this feature. 
In order to achieve a more efficient electrical expression 
of spin accumulation one has to move beyond a passive 
two-terminal device such as spin valve, and consider a 
spin-transistor system in which additional external stim- 
uli can control the magnetoresistive effects. Below we 
review several proposals of devices consisting of more 
than two ferromagnetic terminals connected to a semi- 
conductor channel. Their common feature is the use of a 
ferromagnetic contact kept close to zero bias fSec. Ill Ct . 



which is used to sense the spin accumulation in the semi- 
conductor beneath it. 



A. Magnetic Contact Transistor 

The Magnetic Contact Transistor (MCT) consists of 
three ferromagnetic contacts deposited on top of the 
paramagnetic channel (see Fig. [3^). We concentrate on 
situation in which most of the current driven by the volt- 
age Vl passes between the left (L) and middle (M) con- 
tacts, but other arrangements are possiblci^ The volt- 
age Vu is adjusted to keep the R junction at low biaSf^ 
and we describe its spin-dependent conductance using 
Eq. ([8]). Alternatively,— the R terminal can be connected 
to a capacitor C, which adjusts the voltage of the R ter- 
minal so that there is no net charge current in the steady 
state. We will refer to P (AP) configurations of the L 
and M magnets as corresponding to the spin accumula- 
tion patterns described before. In this section we keep 
the M magnet fixed, and consider the P and AP align- 
ments of the L magnet with respect to M. 

Let us first consider a situation in which the R terminal 
is inactive. Then the L and M contacts constitute a spin 
valve with voltage Vl applied to it. For P and AP align- 
ments of L and M the II current is practically the same, 
but the spin splitting of layer-averaged electrochemical 
potential in the channel Ag varies between two very dif- 
ferent values. Beneath the injecting and extracting con- 
tacts we have \A^^^ /A^^\(x{2L/lf, with the effective 
length of the active channel covered by L and M ter- 
minals l~Wd+WL+WM. The spins accumulated beneath 
the M terminal diffuse out to the right, but if w^-^L the 
spin accumulation beneath the R magnet is practically 
the same as beneath the M contact. 

The operational principle of the "static" MCT relies 
on the fact that we actively control the Vji voltagci^ 
In either P or AP alignment we bias the R terminal so 
that Ir=0. If we then flip the L magnetization, a finite 
Ir current will flow. This is a consequence of large spin 
accumulation in AP configuration, and spin selectivity 



of the R junction {AGr^O). 
current is given by 



Using |Ag^|< 



AP 



•'on 



\Ir\ 



agr Ae 



AP 



A 



the 



(14) 



Thus, we have found a way to digitize the MR effect in 
the R contact. Instead of some finite ratio of P and AP 
currents, we can have zero current for one and a finite 
current for the other configuration. Even after taking 
the voltage noise in the system into account, the "on" 
and "off' currents should be easily discernible for real- 
istic parameters of Fe/GaAs systemi^ This digitization 
of the magnetoresistance had been observed recently in 
MnAs/GaAs three-terminal structure.— 

The digitization holds for the MR effect measured in 
the third (R) terminal. The larger currents in the other 
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FIG. 3: (a) The three-terminal Magnetic Contact Transistor 
(MCT). The outer edges if the channel are removed in or- 
der to confine the spin accumulation under and between the 
contacts.— The R contact is connected either to a control- 
lable voltage Vh, or to a capacitor C which maintains zero 
steady state current in this terminal, (b) A schematic pic- 
ture illustrating the principle of the MCT operation. The 
solid lines are the electrochemical potentials (^a) in the chan- 
nel beneath the R contact for AP and P alignments of the 
L and M terminals. The dashed lines show the values of the 
electrochemical potential inside the R contact for which the 
Ir current is quenched for a given alignment of all magneti- 
zations (represented by three arrows). The energy scale and 
A^'^^/A^^ ratio correspond to an MCT with barrier conduc- 
tances G=10'' f2~^cm~'^, dimensions ui_L=TOi\/=TOj{=400 nm, 
Wd=200 nm, /i=100 nm, and Vl=0.1 V. 



two (L and M) have a small relative change when we 
alternate between P and AP configurations. Yet these 
contacts do almost all of the job of injecting and extract- 
ing spin-polarized currents. We can say that we have 
transferred the magneto-resistive effect to the third con- 
tact, where we can tune it by Vr voltage. We have called 
this s'pin transference.^^ An alternative term of "trans- 
ferable magnetoresistance" underlines the connection to 
standard transistors. 

To sketch the derivation of the above effect we write the 
spin dependent electrochemical potentials underneath 
the R terminal as 



(15) 



where ^ is the average value of electrochemical potential. 
If we apply any voltage Vr, the situation inside the chan- 
nel will change in general. However, we are interested in 



Vr such that there is only a small (possibly zero) Ir 
current. Then the spin accumulation determined by the 
larger II injected into the channel remains practically 
unaffected. For such a voltage applied to R we have the 
current density entering the R contact: 



JR 



Gr I 

e 



i) 



agr Ae 



(16) 



where ^J■R—~eVli, Gr is the total conductance of the R 
junction, and AGr=G^ — G^ is its spin selectivity. The 
corresponding spin current density Aj=j^ — j_ is 

A . AGfl Gr A^ 

AjR = [f^R - TT ■ 17) 

e e 2 

We denote by fiQ the electrochemical potential in R mag- 
net that quenches the total current in this contact: 



AGr AC 
Gn 2 



(18) 



Plugging this /io for one L/M alignment into Eq. ([T| 
with A^ corresponding to the other L/M alignment we 
obtain Eq. (fT4|) . Depending on the alignment of the L and 
R magnets, the Vr voltage which corresponds to Ir—0 
takes on four possible values denoted by dashed lines in 
Fig. When the R terminal is connected to the capac- 
itor, we have four possible charge states of the capacitor 
in the steady state. We discuss the transient currents 
driven by L or R magnetization dynamics for this case in 
Sec. IVB] below. 

In the "off" state, when Vr is adjusted so that the net 
/i^=0, there is a non-zero spin current flowing into the R 
terminal: 



AJr" 



P'^P AGl 



Gr 



2e 



Gr 



(19) 



The electrons with opposite spins flow in the opposite 
directions, giving zero charge current, but adding to a 
net flow of spin polarization. A similar effect was pre- 
dicted in a lateral structure with non-magnetic source 
and drain and two ferromagnetic gates, into which the 
current could leak.'^^ A finite spin current entering a fer- 
romagnet can lead to reversal or precession of the mag- 
netization due to spin-transfer torque if the polarization 
of injected spins is non-coUinear with the magnetization 
axis of the magnet This effect has been observed in 
an all-metallic system^^^ in which the magnetization of 
a floating terminal (zero net charge current through it) 
was switched by a pure spin current. However, in our 
case the relatively low carrier density in the semiconduc- 
tor together with the resistive contacts cannot transfer 
enough angular momentum for such switching to occur. 

A system similar to the MCT has been known for some 
time in the field of all-metal magnetoelectronics as a non- 
local spin valve.— i^^i^ However, the third contact in the 
non-local spin valve is used as a passive floating termi- 
nal (essentially a spin dependent voltage probe). In the 
MCT all the contacts are active terminals controlled by 
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applied voltages. The possibility of control is closely re- 
lated to the use of the non-degenerate semiconductor as a 
paramagnetic channel. Due to very small concentration 
of carriers compared to metals, spin injection can lead to 
spin splittings of electrochemical potentials of the order 
of milivolts in the channel. Then, voltages supplied with 
mV accuracy can efficiently tune the magnetoresitive ef- 
fect measured in one of the terminals. 



B. Electric readout of magnetization dynamics 

The MCT can also be used for electric measurement of 
magnetization dynamics and dynamical readout of mag- 
netization alignment. With the R contact connected to a 
capacitor C there is zero current in the steady state J*^ 
For any alignment of all the magnetizations, the charge 
on the capacitor adjusts itself so that the electrochem- 
ical potential in the R contact is equal to the current- 
quenching /^Oi and there is no need for external voltage 
tuning. When either the L or the R magnet starts to re- 
verse due to an application of an external magnetic field 
pulse, the potential in the R contact changes between two 
steady-state values (dashed lines in Fig. [S]) correspond- 
ing to the initial and final alignments of the magnetiza- 
tions. The measurement of the accompanying transient 
current recharging the capacitor allows for electri- 
cal monitoring of magnetization dynamics. Alternative 
application is a dynamical readout of the L/M alignment. 
In the P case, the 2-k rotation of the R magnet results in 
the lR{t) current oscillation of much smaller amplitude 
than for the AP case (see Fig. |4]) . 

The on-chip manipulation of L and R magnetizations is 
possible using the architecture^^i^l of Magnetic Random 
Access Memory (MRAM). In MRAM the nanomagnets 
(with typical size similar to what we use in our mod- 
eling) arranged in a square array are addressed using a 
network of current-carrying wires positioned above and 
below the magnets. Pulses of current generate time- 
dependent magnetic fields through Ampere's law, and 
these field can be used to switch each of the magnets 
separately. 

The capacitor RC time (with R being of the order 
of the junction resistance, as they are the most resis- 
tive elements in the circuit) has to be at most com- 
parable to the magnetization dynamics time-scale. For 
G'^IO^ ri~^cm~^ and junction area of one /im^, the ca- 
pacitance C=40 fF used below results in RC time of 
about a nanosecond. In all-metallic systems, in which 
the junctions are much less resistive, the RC time is not 
a problem. However, typical spin accumulation in para- 
magnetic metal corresponds to A/i (and consequently 
the voltage swing on the capacitor) of less than 
Then, for the transient current to be measurable one has 
to use a nearly macroscopic capacitor, which rules out 
an application in integrated circuits. Again, the small 
carrier density in a semiconductor (allowing for A/i^lO 
mV) is indispensable. 




FIG. 4: (a) R current signal for reversal of L magnetization 
occurring on time-scale of 3, 5 and 10 ns starting from AP 
alignment of L relative to M magnet, (b) R current signal 
for 2ti rotation of R magnet for P and AP alignments of L 
and M magnets. The period of rotation is 3 ns. The con- 
ductance of the barriers G=10* Q}cm~^ . The area of a junc- 
tion is 1 /im^, and the barrier thickness is taken to be 10 
nm, resulting in junction resisitance i?s = 10 kfi and capaci- 
tance Cb = 10 fP. The external capacitance is C=40 fP. The 
channel is GaAs at room temperature, with carrier density 
no=10^''cm"^. Adapted from Ref. 69. 



In order to model the transient behavior we add the 
time dependence to the formalism of lateral spin diffu- 
sion. We are interested in time-scales of a least tens 
of picoseconds. The fastest magnetization reversal time 
is about 100 ps^^ and magnetization switching times 
used in commercial devices are of the order of a nanosec- 
ond. Thus, we use an adiabatic approximation with re- 
spect to the processes occurring on a much shorter (sub- 
picosecond) time scales: the momentum scattering and 
dielectric relaxation. 

The time-dependent diffusion equation for spin split- 
ting of the layer-averaged electrochemical potential 



9A^ 

dt 



D 



AC 



(20) 



The a(/3) dimensionless parameter is given by 
2L'^{G++{—)G-)/ah. The dynamics of magnetiza- 
tion is parametrized by p{t)^AG{t), which characterizes 
the contact polarization only along the z axis. If we 
deal with a coherent precession of magnetization then 
this is an approximation. In principle, one should treat 
the diffusion of spin accumulation treated as a vector 
quantity;^ and take into account the noncoUinearity of 
spins and the magnets in the tunneling process i^^i^i^ 
However, for tunneling barriers the non-trivial effects of 
this noncoUinearity are expected to be small ^^i^ and 
the only thing that matters is the average polarization 
along the z direction. Then we can model the influence 
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of the contact with magnetization making an angle 6 
with the z axis by assuming that f3^ AG cos 9. On the 
other hand, if the magnetization reversal is incoherent 
(e.g. proceeding by nucleation of domains with opposite 
magnetization), the parameter f3(t) describes an area av- 
erage of spin-selectivity of magnetically inhomogeneous 
contact, and it is proportional to the z component of the 
contact's magnetization. 

The dielectric relaxation (about 100 fs for non- 
degenerate semiconductor with nQ=10^^ cm~'^) is much 
faster than the time-scale of magnetization dynamics and 
spin diffusion, so we assume quasi-neutrality in the chan- 
nel at all times {6n+{t)+6n-{t)—0). In the linear regime 
under consideration (when A^<fcsT) the average elec- 
trochemical potential (, = {£,+ +(,-)/ 2 is equal to —ecj). 
At every moment of time ^ fulfills the Laplace equation 
with boundary conditions given by currents at the inter- 
faces. In the time-dependent case these include also dis- 
placement currents connected with charging of the bar- 
rier capacitance Cb- A Schottky barrier is a dipole layer, 
and its capacitance can have a strong effect on dynam- 
ics of currents on time scales of interest here. With the 
displacement current taken into account, the boundary 
condition for spin current is: 

J, = ^int) ~ m) + ^^i/^^w - ' (21) 

where cb is the barrier capacitance per unit area. The 
second term in the above equation represents the carriers 
which flow towards the barrier, but do not tunnel through 
it. Instead, they stay in the semiconductor close to the 
barrier, making the depletion region slightly thinner or 
wider. The charge involved in this process is negligible 
compared to the charge already swept out from the semi- 
conductor, so we can keep cb constant. For small spin 
splitting (so that the conductivities cr+~cr_) the same 
amount of carriers of each spin is going to be brought 
from the channel into the barrier, and the displacement 
current is the same for each spin in Eq. (|2ip. For layer- 
averaged ^ we get then 

with the right hand side of Eq. ([^^ non-zero only under 
the contacts. The magnetization dynamics of i*^ con- 
tact translates into time-dependence of /3i, driving the 
spin diffusion in Eq. (pOI and electric potential in the 
channel in Eq. (P^ . From we calculate the current 
//{(t) charging the capacitor C. The electrochemical po- 
tential of the R terminal Mi?,— "^Vr changes according to 
dVul dt—Ifi/ C. Examples of calculations for two possible 
modes of operation (sensing the L dynamics and reading 
out the L/M alignment) are shown in Fig. U) 

C. Reprogrammable Magneto-Logic Gate 

The same physical principle of operation can be har- 
nessed to achieve a higher level functionality. In Fig. [5^ 





a 5 10 15 

time (ns) 



FIG. 5: (a) A five-terminal magneto-logic gate (MLG). The 
logic inputs '0' and '1' are encoded by magnetization direction 
of the A,B,X,Y terminals (see text for details). As shown here, 
the gate is set to work as a NAND operation between X and 
Y (A and B are fixed to '1' values). In the read-out phase the 
magnetization of the middle (M) terminal is rotated by 2-k, or 
the back-gate voltage Vg is increased, (b) The lM{t) transient 
current triggered by the rotation of the M magnetization. The 
small signal for X='l', Y='l' corresponds to logical '0' output. 
Adapted from Ref. 38. 



we present a scheme of a five-terminal system,~ in 
which the electric sensing of spin accumulation is used 
to perform a logic operation, i.e., two bits of input 
are converted into a binary output signal. This is a 
reprogrammable magneto- logic gate (MLG). Spintron- 
ics logic gates have been proposed in purely metallic 
systemS j^^i^^'^'^i^^'^^i^'^ but ours is the first proposal which 
employs semiconductors as active elements of the system. 

The system presented in Fig. [5^ works in the following 
way. The charge currents are flowing between two pairs 
of terminals (X and A, Y and B), between which the bias 
Vdd is applied. Depending on the alignment of these pairs 
of magnets, different patterns of spin accumulation are 
created in the channel: if both X/A and Y/B are AP the 
spin accumulation underneath M is large; if only one pair 
of contacts is AP the spin accumulation is approximately 
two times smaller, and if both pairs are P there is a very 
small beneath the M terminal. The M contact is 
used to directly express the differences in the average 
spin accumulation beneath it. 

The logical inputs are encoded by magnetization direc- 
tions of A, B, X, and Y terminals. We will concentrate 
on the case in which A and B magnetizations are pre- 
set, defining the logic function of the gate. This repro- 
grammability is an important feature of magnetization- 
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based logic. X and Y are then the logic inputs, and 
the output is generated when the M magnet is rotated 
by 27r, triggering a transient /^/(i) current of amplitude 
proportional to the spin accumulation in the middle of 
the channel. Let us focus on the example of the NAND 
gate, as any other logic function can be realized by using 
a finite number of such gates. For the NAND opera- 
tion, A and B magnets are set parallel to each other in 
the direction defining the logical '1'. The amplitude of 
the lM{t) oscillation is two times larger for ^='0' and 
y='0' compared to the case when one of them is '0' and 
the other is '1', and for '11' case the current is negligi- 
ble. This is shown in Fig. [SJd. The transient current can 
be captured by an external electronic circuit, and then 
used to control a suitable write operation applied to the 
magnetic contact of another gate.^^ 

Instead of using magnetic field pulses to drive the 2ti 
rotation of the M magnet, one can employ the idea of the 
spin switch outlined in Sec. Ill Bl in which the M magne- 
tization is pinned, and the profile of the conduction band 
beneath M is changed by applying a voltage pulse Vq (t) 
to the back gate shown in Fig. [5^. Vg is chosen to deplete 
the electrons from the lower part of the channel, biasing 
them vertically towards the Schottky barrier. Thus, the 
effective spin selectivity of the M junction is switched as 
discussed in Sec. IIIBI This is qualitatively the same as 
reversing the M magnetization, and the transient cur- 
rent Im should be generated. However, its quantitative 
calculation is much more involved than in the case of 
magnetization rotation, and it will be addressed in fu- 



ture work. 



V. SUMMARY 

We have reviewed the physics of spin injection and ex- 
traction through inhomogeneously doped Schottky barri- 
ers. We have paid special attention to the states localized 
close to the junction as a consequence of heavy doping 
near the interface, and we have shown that they play a 
crucial role in the recent spin extraction experiments. A 
junction held at low bias in the presence of spin accu- 
mulation in the semiconductor has been also analyzed, 
and how its spin-dependent properties could be used to 
electrically sense the spin accumulation in a semiconduc- 
tor has been described. The multi-terminal systems (a 
three-terminal transistor and a five-terminal logic gate) 
which we have described rely only on spin-selectivity of 
the junctions and the presence of the spin accumulation. 
As such, they are ideal candidates for spintronics devices 
working at room temperature and, possibly, in silicon. 
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